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We use renormalization group to calculate the reunion and survival exponents of a set of random 
walkers interacting with a long range 1/r and a short range interaction. These exponents are 
used to study the binding-unbinding transition of polymers and the behavior of several quantum 
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I. INTRODUCTION 



The problem of reunion of interacting random walkers 
appears in many situations, directly or in disguise, such 
as e.g., in situations involving one-dimensional string- 
like objects such as interfaces, steps on surfaces, in one- 
dimensional quantum problems, in several types of poly- 
mer problems, wetting, etc. The issue quite often is the 
large-length scaling of the partition function of reunion 
(or survival) of a set of walkers all starting, say, at origin. 
For a purely short range repulsive interaction, this is the 
problem of "vicious walkers." Many of the vicious-walker 
reunion and survival exponents are known from exact 
calculations in one dimension |j]J^], the renormalization- 
group (RG) approach for general dimensions, and 
lattice models in two dimensions j5[. Recent studies of 
phase transitions on vicinal or miscut surfaces re- 
quired the scaling behavior of such partition functions in 
the presence of a long range interaction (namely 1/r 2 ) 
in addition to short range interactions. The effective in- 
teraction of steps on a vicinal surface has a long r~ 2 tail 
H and recent studies on Si(113) surfaces indicate the 
presence of an attractive short range interaction || as 
well. Similar long range interactions occur in the one- 
dimensional Calogero-Sutherland model of interact- 
ing quantum particles, and as the angular momentum 
term in higher-dimensional quantum problems. The re- 
union and survival behaviors turn out to be the unifying 
feature in these wide varieties of problems. These moti- 
vated us to study the scaling limit of various partition 
functions in the presence of long range interactions, es- 
pecially r~ 2 . We use these results for studying phase 
transitions of polymers and then show how many results 
of quantum problems can be recovered from such an anal- 
ysis. 

The success of the RG approach for short range in- 
teractions P-fslJrH] and the possibility of exact renormal- 
ization in the presence of long range interactions Jl2||7| 
allowed us to use the renormalization-group approach for 



this problem. It is known that the probability distribu- 
tion function for a random walker can be represented 
by the Wiener integral, which as a path integral could, 
a la Edwards, represent a polymer. Any interaction of 
two or more walkers at the same time or step length 
would translate into an interaction among the polymers 
at the same contour length. For d-dimensional walkers, 
we would then get (d + l)-dimensional interacting di- 
rected polymers. We shall be using this polymer language 
throughout this paper. 

Phase transitions in polymerlike systems by interac- 
tions at the same contour length occur, e.g., in the case of 
steps already mentioned, in flux lines in superconductors 
Jl3f , in melting and unzipping of DNA |i4-|l9|], to cite a 
few. A directed polymer formulation then becomes quite 
natural. In this paper, our focus is mainly on the effect 
of the long range interaction on these phase transitions, 
recovering the results for the short range case as a special 
case. Furthermore, we exploit the relation between poly- 
mers and quantum problems to show the importance of 
reunion behaviors in wide varieties of problems in various 
dimensions [ p0[ . 

In this paper, we recover several results known from 
exact solutions, illustrating how the reunion exponent 
armed with the renormalization group gives a unified 
approach to these problems. The RG approach can be 
extended to other types of interactions, especially other 
marginal interactions, for which exact solutions are not 
known, such as that required in the vicinal surface prob- 
lem of Ref. @. There lies the merit of the approach of 
this paper. 

The Hamiltonian and the general exponents of inter- 
est are introduced in Sec. II. The general RG approach 
is discussed in Sec. Ill, where the exponents are also cal- 
culated, with some of the details in the Appendix. These 
results are then used for polymer problems in Sec. IV 
and to quantum problems in Sec. V. The conclusion is 
given in Sec. VI. 
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II. HAMILTONIAN AND THE EXPONENTS 

The Hamiltonian for p-directed polymers with pairwise 
interaction MM is 



(2.1a) 



where r is a d dimensional position vector and 
r,(z) — Yj(z). The interaction potential is taken as 



V(v) = v 6 A (r) + 



(2.1b) 



subsequent sections on polymers and quantum problems, 
we need mainly We therefore adopt the notation 



* = * 



R,2- 



(2.3) 



Our primary interest is in the marginal case a = 2, 
though other values are also discussed briefly. The im- 
portance of er = 2 can be gauged if the Hamiltonian is 
considered as a quantum-mechanical system in imaginary 
time t — iz. The kinetic energy and the r~ 2 interaction 
(like the angular momentum) scale in the same way for all 
dimensions. The r -2 potential is the centrifugal barrier 
in quantum mechanics, while it occurs in surface prob- 
lems as an effective interaction between steps induced by 
the elastic strains on the surface Hi. 



and <5a (r) is a delta function with a cutoff such that in 
Fourier space 



<5 A (q) = 1 for q < A, 
= otherwise. 



(2.1c) 



We shall use this cutoff for the whole potential. Our 
interest is in the asymptotic behavior of the following 
partition functions. 

(a) Reunion of all the chains of length N at a point r, 

= J VR e- H - f[ [6 d (r t (0)) S d (r t (N) - r)] (2.2a) 

1— 1 

N-^"; (2.2b) 

(b) reunion anywhere, 

(N) = J d d r Z R , p (N, r) (2.2c) 

^ n -*r,p. ( 2 . 2d) 

and (c) survival, 

Z S , P (N) = JvR e~ H * n[^(r,;(0))] (2.2e) 



pj--4>s, P 



(2.2f) 



where J VR stands for the summation over all possible 
configurations of the polymers (path integrals). The con- 
straint that all chains are tied together at the origin is 
represented by the product of the 5 functions. In case 
(c), the ends at z = N are free for all the chains. 

In the above equations, the asymptotic behaviors de- 
fine the basic exponents (ipR^,^ R, p ,ips,p) of interest in 
this paper. These exponents are expected to be univer- 
sal, independent of the detailed microscopic form of the 
polymers. Hence the choice of a continuum model. 

The notation introduced above is used in Sec. Ill, 
where these exponents are calculated. However, in the 



A. Known exponents for the short range case 

For the noninteracting problem, vq = g = 0, the ex- 
ponents follow from the Gaussian behavior of the chains. 
The partition function for a Gaussian chain of length N 
is given by 

Z(RiV|00) = (2itN)- d/2 exp(-i? 2 /2iV), (2.4) 
from which one obtains 



-0S,p = 0, IpR, 



pd 



and * 



R./j 



(P - l)d 



(2.5) 



It is also known that for purely repulsive walkers (vicious 
walkers) in d = 1 (i.e., with g = 0) 

VS, P = 7 , VR, P = y, and * R:P = — - — . 

(2.6) 

We earlier generalized the above results to d = 2 — e 
for the short range case (i.e. g — 0), by using RG and 
obtained ||,| 

. , pd (p— l)d 
VS,p = Vp, WR, P = y + 2^, and W RjP = h 2r) p , 

(2.7a) 



where, with e = 2 — d 
2rj p 



2> + 3 3 ln(3 ^ + 



(2.7b) 



For d = 2, our RG yields Z RiP - ^-^-^(lnTV)^^" 1 ), 
which has since been obtained by the exact lattice calcu- 
lations of Essam and Guttmann |i| for p = 2. This lends 
further support to the validity of the renormalization- 
group approach. 

That only one anomalous exponent is needed was 
shown explicitly in Refs. This could be understood 

from the geometry that a reunion partition function can 
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be thought of as two survival partition functions con- 
nected together. 

In this paper, we calculate the anomalous part rj p for 
the case with long range interaction at the one loop level. 
Since this one- loop RG generates exact results J12], we 
obtain exact exponents with nonzero g. 

III. RENORMALIZATION GROUP 



u = vL c , g = K d Agl? a , and u = K d u H — — - — (3.5) 

a — a 

with v — vo as the bare value of the short range coupling 
constant. Note that for a = 2, g = g. 

By implementing a thin-shell integration over A—dA < 
k < A and then rescaling the cutoff back to A (= 1), the 
coupling constants' renormalization can be written in the 
form of recursion relations as 



A. Survival partition function 

To evaluate the survival exponent, we use a 
momentum-shell RG approach as used in Ref . (l2) . The 
renormalization process requires integrating out small- 
scale fluctuations, so that the effect of interactions of the 
two chains within a small distance is taken into account 
by redefining the parameters of the Hamiltonian. In this 
problem, in addition to the parameters of the Hamilto- 
nian, we need to consider the survival partition function 
also. 

In Fourier space, the renormalization of the interaction 
can be carried out using the effective interaction at the 
one-loop level 



V eS (k) = V(k) 



d d q 7(|k-q|)V(g) 



(2tt) c 



The Fourier transform defined generically as 



F(k) 



d d r e ik r G(r) 



(3.1) 



(3.2) 



requires an appropriate analytic continuation of the in- 
tegral for singular potentials. In particular, the Fourier 
transform of r~ CT has a singular part 



A — in general, 

a — a 



where 



A = 2 d - a+1 TT d/ 



2 T((d- a + 2)/2) 



(3.3a) 



(3.3b) 



T(a/2) ' 

and, for a = 2, this singular part takes a simpler form, 
1 k 2 - d 



K d d-2' 



(3.3c) 



where K d = S d /(2%) d , S d = 2vr d / 2 /r(d/2) being the sur- 
face of a (i-dimensional unit sphere. A natural choice for 
the potential is therefore 



V(k) =v + 



A 



d — a 



g k 



a—d 



(3.4) 



Let us also introduce at this point the dimensionless pa- 
rameters 



du 2 
L— =eu-u + g, 
dL 



(3.6a) 
(3.6b) 



The thin-shell integration is analytic and therefore does 
not renormalize a singular g-type long range poten tial in 
the Hamiltonian. The flow equation of Eq. ( 3.6b ) for g 
then follow s from dimensional analysis. The exactness of 
Eq. ( 3.6a ) for g = and for g ^ is discussed in Refs. 
[UlLpJ- This marginality of g is a reflection of the impor- 
tance of an r~ 2 potential in the Hamiltonian, as is well 
known in quantum mechanics. In the surface context, 
this marginality implies that the strength of interaction 
between two steps is invariant under scale transforma- 
tion. Henceforth we consider a — 2, and use g instead of 
9- 

The flow equation for u gives rise to two fixed points 
1 



<„ = -[ e ± v / ^+43], 



(3.7) 



where the subscripts s and u represent stable and unsta- 
ble fixed points, respectively. 





(a) 



(b) 



mukherji-bhattacharjee 
FIG. 1. Zeroth- and first-order diagrams for survival. The 
solid lines indicate the polymers and the wavy line interaction. 
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For the survival partition function, the one-loop con- 
tribution is shown in Fig. [l[ The loop contribution is 
similar to the one-loop case for the interaction, though 
here it is 0(u). The flow equation for Zg is (see Appendix 
A for details) 



T dZ s 
~dL = ^ 



so that at the stable fixed point of Eq.(3.6a) 



m = u*/2. 



(3.8) 



(3.9) 



* = * 



R,2 



d 
2 



+ V e2 + 4 .9 



= 1 



\A 2 + 4 .9 



(3-12) 



This is used in the following sections. 

A few things are to be noted here. First, the fixed 
points become complex for g < — e 2 /4, in particular 
for g < —1/4 for d = 1. For this borderline case 
d = 2 ± 2 y/—g, the reunion exponent is = 1. The 
locus of the points in the d vs g plane for $!r 2 = 1 an d 
^R,2 = 2 are shown in Fig. ^. As we show in Sec. IV, 
these lines are very special. 



This result can be generalized to any number of chains 
for which there will be a combinatorial factor m?j from 
the pairwise interaction to yield 



Vp 



Ifp 



2 \2 



(3.10) 



This exact relation for p = 2 is a consequence of the 
connection between the coupling constant u (vertex func- 
tion) and the survival partition function. The interaction 
u can be thought of as two p = 2 survival partition func- 
tions joined at the node so that renormalization of u is a 
product of the renormalizations of two survival partition 
functions. 



first order 




-l u - 3 g 

FIG. 2. The locus of points in the d vs g plane for re- 
union exponent ^ r,2 = 1 and ^ r,2 = 2. The shaded re- 
gion is the region for continuous transition. The transition 
is first-order outside it. The lightly shaded region has com- 
plex-conjugate fixed points that are of importance in the con- 
text of non-hermitian Hamiltonians. 



IV. BINDING-UNBINDING TRANSITION OF 
POLYMERS 

Let us cons ider two directed polymers, the p = 2 case 
of Eq. (2.1a). The fixed-point diagram of Fig. IVA 
shows that at d = 1, for a given g > —1/4, there will 
be a binding-unbinding transition as v is changed (see 
Fig. 1 of Ref. @ ) . The transition is defined by the un- 
stable fixed point while the stable fixed point describes 
the "high-temperature" phase. Several features of this 
transition have been discussed in Ref. |fl2fl in the renor- 
malization group framework and in Ref. |21|], We show 
here that certain unresolved issues (e.g., order of transi- 
tion) in the RG framework can be sorted out by using 
the reunion exponent, thereby gaining a complete pic- 
ture of the problem. The problem can be anticipated by 
noting that the exponent v associated with the diverg- 
ing length scale, as the transition is reached, approaches 
the limit v = 1/2 as g — > 3/4. Care should, therefore, 
be exercised in drawing conclusions from the flow dia- 
gram. We discuss now how the behavior at the stable 
fixed points determines the phase transition at the un- 
stable fixed point. 

For generality, we develop the procedure in a general 
way and then derive the results for the 1/r 2 interaction, 
especially for p = 2 (two-chain case). The internal con- 
sistency is also shown by a finite-size scaling argument. 
That the exact results are recovered corroborates the util- 
ity of the RG approach, coupled with the reunion behav- 
ior, for other interactions as well. 



The reunion-anywhere exponent for d = 1 is given by 



A. General approach 



p-l 



(3.H) 



A special case of this is *f> = \ + u*. For only short 
range interactions (i.e., g = 0), u* = 1 and one gets 
back the vicious walker exponent of Huse and Fisher 
[\& Rp = (p 2 - l)/2] Ql]. One notes that the anoma- 
lous part is just the combinatorial factor. The general 
results for all three reunion and survival exponents are 
obtained by substituting Eq. ( 3.10j ) for rj p in Eq. (2.7a). 
We quote the exact result for d = 2 — e and p = 2, 



In order to study the nature of the transition, we adopt 
the procedure of Ref. JlJ . The partition function is ob- 
tained by summing over all configurations of the two 
chains. Our main interest is in the phase-transition be- 
havior and therefore only long-distance behaviors will be 
considered. The configurations can be characterized by 
the alternate sequence of domains of bound (A) a nd un- 
bound (bubble, B) regions as shown Fig. IVA. The 



bound regions correspond to the left side of the unsta- 
ble fixed point of the flow diagram (the coupling con- 
stant flows to — oo), while the bubbles correspond to the 
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high-temperature phase characterized by the stable fixed 
point. 




(b) 



transition B 
point 

• — > 



u u u s 



Mukherji-bhiittacharjee. Fig .1 

FIG. 3. (a) A general configuration of the two polymers, 
consisting of bound (A) and unbound (B) regions, (b) A 
generic fixed-point diagram for u from the RG flow. See also 
Fig. f of Ref. [7]. 



Let us denote the partition function for a segment of 
length N of type X by Q^. We assign a weight v to each 
A-B and B-A junction, and for simplicity we consider 
cases in which the chains start and end in A-type do- 
mains. This extra factor v is to take care of the deviation 
from the strict A- or B-type behavior at the junctions. 
The total partition function can be written as 



Qn ~ Q 



N 



dz2 / dz\Q 



N-Z 2 l 



+ 



(4.1) 



To exploit the convolution nature of the terms of the se- 
ries, we go over to the grand-canonical ensemble with 
respect to the chainlength. This is equivalent to taking 
the Laplace transform of the partition functions with re- 
spect to N. Defining the Laplace transform as 



G x ( s ) 



e~ sN Q% dN, 



(4.2) 



where X is A, B or T, the above series can be written 
in the form of a geometric series and can be summed to 
obtain 



G T (s) 



l-v 2 G A {s)G B {s)' 



(4.3) 



The existence of the thermodynamic limit ensures 
that G x (s) is analytic for s > fx, where fx — 
limTv^oo InQiV) an d G x (s) — > as s — > oo. In 
other words, the largest real singularity of G x (s) in the 
complex-s plane determines the corresponding limiting 
free energy per unit length p2fl . For G T (s), this could 
either be the singularity of G or G B or may come from 
the vanishing of the denominator of Eq. ( |4.3| ), i.e. the 
solution of 



v 2 G B {s) = l/G A (s). 



(4.4) 



I n ca se G A or G B diverges, then the denominator of Eq. 
(4.3) vanishes for a larger value of s, and then it is the 
root of Eq. (4.4) that is the relevant one. In such a 
situation, the free energy will not have any singularity. 
It follows that it is only the nondiverging branch points 
that could lead to phase transitions. 

For the A-type region, the partition function = 
e -f3eN ^ w j iere p ig inverse temperature. Its Laplace 
transform has a simple pole. Therefore, the singularity of 
G A is not important for our discussion. Consequently, it 
is the partition function of the bubbles (high-temperature 
phase) that determines the behavior of the chains [^3] . 
The partition function of the bubble at the stable fixed 
point for large N is given by 



Qn p ' 



-Na {T) 



AT* 



(4.5) 



where the reunion-anywhere exponent at the stable 
fixed point appears. Here and in the following discus- 
sion, for convenience, we have removed the subscripts 
of th e reunion exponent introduced in Sec. II [see Eq. 
fl2.3| )] . The Laplace transform of this partition function 
has a singularity at s = ctq (T) , and the nature of the sin- 
gularity depends on the value of the reunion exponent, 
"J". There is a divergence at this singularity if "J < 1. 
A phase transition therefore occurs only for \& > 1. The 
details of the graphical solution, for the discrete case, can 
be found in Ref. 0. 

For \P > 1, the h igh temperature phase is described 
by the root of Eq. (4.4) and the transition is give n by 
the temperature T = T c at which the root of Eq. (4.4) 
coincides with the singularity of G B (s) . Defining the de- 
viation of temperature as t = T c — T, the free energy 
[from the root of Eq. 4.4] can be written as 



/t w 0o (T)— I t [VC*-!) for 1< * < 2, (4.6a) 

rn— 1 

« a (T) - £ aj t?+ | t 



for m < ^> < m + 1 . 



(4.6b) 



We find a critical behavior for 1 < \& < 2 but a first- 
order transition for ^ > 2. In case of a critical behavior, 
the specific-heat exponent can be read off from the free 
energy as 



a = 2- l)" 1 = 



2* 



* - 1 



(4.7) 



B. d=l 

The crucial feature in the approach developed in the 
preceding subsection is the alternate sequence of the two 
types of "bubbles" and therefore the results can be used, 
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for example, in the case of unbinding of p chains with 
p-body interaction and other cases. To take advan- 

tage of previous studies, we concentrate on the two-chain 
problem. 

Let us now consider the case of two directed poly- 
mers in (1 + 1) dimensions, for which we ob tained \E' 
i I + y/l + 4g)/2. The fixed points (see Fig. [IV Aj ) indi- 
cate that there is a binding-unbinding transition at the 
unstable fixed point for any given g as u is varied. For 
any uq < «*, the renormalized flow goes to — oo, so that 
a length scale can be identified from u(l*) — > — oo. From 
this, as u — > u*, the diverging length scale comes out as 
(t = Au = u — u*) 



For c? = 1, we have 



with u± 



1 



1 



(4.8) 



(4.9) 



This shows that v± = 1/2 at g = 3/4, a signature of 
something special happening there. The reunion expo- 
nent, coupled with the analysis of the preceding subsec- 
tion, tells us that there cannot be a phase transition for 
g < —1/4, a fact corroborated by the fixed-point diagram 
(see Fig. 1 of Ref . [f? | ) . A critical behavior is expected for 
3/4>.g>— 1/4 and a first-order transition for g > 3/4. 

For —1/4 < g < 3/4, which includes the pure short 
range g = case, the average fraction of the length in 
the bound state is given by 



0(t) 



df 
dt 



t 



(4.10) 



This defines the order-parameter exponent j3. This ex- 
ponent and the specific-heat exponent are given by 



* - 1 



and a = 2 



-1 



(4.11) 



Note that a = 1 at g = 3/4, a requirement for a first- 
order transition. 

The bubble lengths have fluctuations and this gives the 
measure of the diverging length scale parallel to the chain 
as 



en ~i*r 



with 



2i/i 



(4.12) 



Note that the hyperscaling relation dv\\ — 2 — a is obeyed 
by these nonuniversal exponents with d = 1. This is 
because the free-energy density is the free energy per 
unit length of the polymers and not the unit (d + 1)- 
dimensional area. 

If g > 3/4, 0(t) — > const, as t — * 0—. This indicates 
a first-order transition. The longitudinal length-scale ex- 
ponent is z/ii = 1, i.e., it sticks to its value at g = 3/4. 



However, the free energy as given by Eq. (4.6a) shows a 
weak singularity that will be reflected in the divergence 
of an appropriate higher derivative. One therefore finds a 
rather unusual first-order transition with weak singular- 
ities and diverging length scales pl[ . Such scales can be 
determined from the higher cumulants of the length fluc- 
tuations of the bubbles. We discuss these issues in the 
next section in the context of the equivalent quantum 
problem. 

^ increases with g and every-time \P crosses an integer, 
the diverging derivative shifts by 1. These special values 
are g = (m — l) 2 — 1/4. 

At g — —1/4, two fixed points u* and u* merge at a 
common value u* — 1/2. The bound phase corresponds 
to the uq < 1/2 region and the transverse length scale 
diverges as 



a~exp[l/(l/2-tio)]. 



(4.13) 



However a Kosterlitz-Thouless-type behavior is observed 
if g is varied fl^| . 



1. Finite size scaling for the critical point 



We have derived the exponents for the critical point 
from our results of the stable fixed point or the high- 
temperature phase. An independent verification of the 
exponents comes from a finite-size scaling argument used 
strictly in the critical region characterized by the unsta- 
ble fixed point. The reunion describes the contacts of 
the two chains and therefore the number of contacts at 
the unstable fixed point would have a finite-size scaling 
behavior (L is a length along the chain) 



e. 



which identifies 



(4.14) 



2 - vfM ij/M 

V]/{s} _ i " 



'II 



(4.15) 



where we have introduced the superscripts {s} and {u} 
to distinguish the values of the exponent at the stable 
and unstable fixed points, respectively. of Eq. ( 4.11 ) 
isjFjj 5 } here.] Using the unstable fixed point value in Eq. 
fl3.ll ), we do see this equality to be true because 



(4.16) 
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For 1 < d < 2, i.e., d — 2 — e with positive e, the results 
of the p reced ing subsection can be repeated by using 
of Eq. ( 3.12 ). The reunion exponent requires g > — e/4 
(so that W > 1) for a phase transition that is identical 
to the condition of the real fixed point of the flow equa- 
tion. The length-scale exponent becomes equal to 1/2 at 
g = d(4 — d)/2 = (4 — e 2 )/2, which coincides with the 
value of g at which \& = 4. Th e exp onent s sat isfy the 
general scaling relations of Eqs. (4.12) and (4.15). 

The stability of the fixed points flips at d = 2. The 
results for d > 2 can be obtained by doing an analytic 
continuation from e > to e < 0. For example, on the 
stable branch now '5 = 1 + \J e 2 + Ag/2. Figure]^ shows 
the regions of critical and first-order transition in the d- 
vs-g plane. 



V. QUANTUM PROBLEM 



The Hamiltonian of Eq. (2.1a) represents two quan- 
tum particles either in the path-integral representation 
in imaginary time or as a statistical mechanical problem 
at T = 0. In the center-of-mass frame, the equivalent 
Schrodinger equation for a = 2, where vq(t) is the short 
range part of the interaction, is 



V 2 ^(r)+ u (r) + 4 ^(r)=£ty(r). 



(5.1) 



For zero energy, the radial part of the wave function can 
be written as 



R{r) = exp 



u(r) 



dr 



with u(r) satisfying 



d)u - 



(5.2) 



(5.3) 



Taking I = ln(r/a), we can recast the above equation in 
a form resembling the RG flow equation. In other words, 
the RG flow equation at long distances determines the 
radial part of the zero-energy wave function. 

We consider the transition point itself, which corre- 
sponds to the unstable fixed point. In case u of Eq. ( |5.3| ) 
reaches a fixed point, the wave function has an algebraic 
tail as R{r) ~| r 

The wave function R(r) satisfying the Schrodinger 
equation is analogous to the partition function of a poly- 
mer starting from origin to the point r. Its behavior 
is similar to the survival partition function because in 
the process the particle may feel (in a perturbative ap- 
proach) the potential any number of times. The scal- 
ing that time and space are related by [N] = [r] 2 then 
gives, via the scaling of the survival partition function, 



1o r 



|2*s 



This is what we get from Eq. 



At d = 1, we th en get at u = u* [the unstable fixed 
point of Eq. (|/7j) with e = 1] R(r) ~| r \<. The mo- 
ments of scalar r 2 = r • r, defined for general d by 



< r 2 P >= 



J(r ■r)PR{r) 2 r d ~ 1 dr 
Ji?( r )2 r rf-i dr ' 



(5.4) 



would depend on the long tail of the wave function. The 
wave-function itself is not normalizable for u* > —1/2, 
i.e., for g < 3/4. This is a signature of an unbound state. 
In such a situation, by tuning the short range part of 
the potential from below, one can get a bound state ar- 
bitrarily close to zero energy. As the zero-energy state 
is reached, the length scale that measures the bounded- 
ness or localization of the state increases without bound. 
Obviously, all the moments are divergent. This is the 
quantum picture of the criticality discussed in the pre- 
ceding section for polymers. 

For g > 3/4, the zero-energy wave function is normal- 
izable. This means a bound state with zero energy but 
with a long tail of the wave function. In this region, the 
bound state can also be made to approach the zero energy 
state but there will always be a finite length scale com- 
ing from the finite moments of the wave-function. This 
is a first-order transition, but with a difference. The mo- 
ments of r 2 are finite for p < (y/l + Ag) /2 — 1, i.e., for 
p < ^ — 2, and any value of p < ^ — 2 may be used 
to define a finite length scale. However, the divergence 
of the higher moments indicates some remnants of the 
"criticality." The integer values of the re unio n exponent, 
which show up in the free energy in Eq. (4.6a), also make 
their presence felt here as the special values at which a 
new integer moment becomes finite. 

We immediately see that all the moments of r 2 are di- 
verging for u* > —1/2, and the first p moments are finite 
for 



g> (p+ l) 2 -l/4. 



(5.5) 



For p violating this criterion, using a finite-size scaling 
analysis (i.e., by cutting off the integral for large r by 
the length scale), one gets 



1 . 



2-^TTIg). 



(5.6) 



That there is a diverging length scale coming from higher 
moments gives an a posteriori justification of the RG 
analysis based on fixed points. 

The states with E < always have at least exponential 
(or faster) decays of the wave function at large distances. 
This rapid decay ensures finiteness of all moments. How- 
ever, at the transition point, the possibility of diverging 
moments arises because of the power-law ( "critical" ) de- 
cay of the wave function. 
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A. Short range interaction, general d 



VI. CONCLUSION 



We now show that the one-dimensional exact results 
for the reunion behavior can be made to bear upon the 
higher-dimensional problems also. In the process, we re- 
cover many results obtained earlier from detailed exact 
solutions for each d by using the properties of special 
functions |§f||f). 

Let us consider a short range central potential V(r). 
Defining <ft{r) — r^~ 1 ^ 2 i?(r), where R(r) is the radial 
part of the wave function, the Schrodinger equation in d 
dimensions can be written as [E6| 



<P±_ t fMd) 

dr 2 



+ V(r) cj>{r) = Efa 



(5.7) 



where Ai(d) = (d + 21 - 3)(d + 21 - l)/4 is the coeffi- 
cient of the angular momentum ( or centrifugal) barrier, 
I being the integer angular momentum quantum number. 
Note that the factor r^ -1 '/ 2 makes the integral over r 
the same as that of the one-dimensional problem. The 
(i-dimensional problem is then reduced to the one dimen- 
sional problem and the exact results for d = 1 for various 
g as obtained in the previous subsection can be used to 
get the features as d and I are changed. The ground state 
is obtained from I = for d > 1. For other values of d, 
I is to be chosen such that 1 — 2l<d<3~2l. (The 
choice makes the minimum centrifugal barrier.) By us- 
ing our results, we see that for Ai(d) < 3/4, there is no 
zero-energy bound state, though the state can be reached 
continuously from below. Substituting the expression for 
Ai(d), we find * = | + I for d > 2. This shows that for 
I = and 1 < d < 4, there can be no bound state at 
zero energy. As the short range parameter is tuned, the 
localization length of the bound state diverges with an 
exponent 



1 



1 



d-2 



for 1 < d < 4, 



(5.8) 



with = 2 at d = 4. For d < 1, the ground state comes 
from I = 1, so that 



1 

~~d 



for 



1< d< 1. 



(5.9) 



In conclusion, we have investigated the reunion of ran- 
dom walkers having both short range and long range in- 
teractions. By using a momentum-shell renormalization- 
group technique, the reunion and survival exponents have 
been calculated. The exponent ^> = ^>r^ for polymers 
has been evaluated at the unbound phase and at the 
binding-unbinding transition point represented by the 
stable and the unstable fixed points, respectively, in the 
coupling constant space. The value of this exponent in 
the unbound phase is crucial in determining the nature 
of the binding-unbinding transition p^,ff2,R|. This tran- 
sition is critical for 1 < 'f < 2 and first order with higher 
moments diverging for > 2. Since ^ is explicitly de- 
pendent on the strength of the long range interaction g, 
the order of thephase transition depends on this param- 
eter. See Fig. g. For example, at d = 1, the dimension 
more relevant in the context of the experimental obser- 
vation of vicinal surfaces Q , one finds a first-order tran- 
sition for g > 3/4 and criticality for —1/4 < g < 3/4. In 
the quantum-mechanical picture, the different nature of 
the phase transition is reflected in the approach of the 
bound state to the zero-energy state as the short range 
part of the potential is tuned. For g < 3/4, one can get 
the bound state arbitrarily close to the zero-energy state 
with a diverging length scale as the gap vanishes, whereas 
for the first-order case (g > 3/4), the length scale remains 
finite. In the latter case, it is possible to defin e di verging 
length scales from higher moments [see Eq. ( |5.6| )]. This 
one-dimensional case can further be extended to the case 
of a quantum particle subjected to a short range potential 
in dimensions d ^ 1 with the centrifugal barrier playing 
the role of the long range potential. The fact that all the 
details and the nuances of the r~ 2 interaction problem 
known from exact solutions Jl0| , pl] , p6| , p7| could be recov- 
ered in a unified manner via the reunion behavior lends 
credence to the general approach developed in this pa- 
per. Our method can be used for other problems and 
interactions as well. 



APPENDIX A: RENORMALIZATION OF THE 
SURVIVAL PARTITION FUNCTION 



Similarly, u±_ = 1/ | d + 2 | for -3 < d < -1 when I = 2. 
In general, 



1 



d + 2l 



for 1 - 21 < d < 3- 21. (5.10) 



For all d < 4, v\\/v±_ = 2. For d > 4, we find from 



Eq. (5.5) that the pth moment of r will be finite 
p < d — 4, and for p > d — 4 one gets 

1 rf-4 

V-n = • 

p 2 2p 



if 



(5.11) 



In the noninteracting case, the survival partition func- 
tion Z s for Gaussian walkers is unity. The anomalous 
exponent of N in the survival partition function appears 
due to the interaction among the chains. The nontriv- 
ial contribution of this interaction is apparent from the 
one-loop 0(u) term, which is shown in Fig [y. The con- 
tribution of this diagram is 

J J dridriZiz, ri | 0, 0)Z(z, r 2 | 0, 0) x 



(v S(ti - r 2 ) 



\ri-r 2 I' 



-) (Al) 
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In the Fourier space, this term appears as 



JV 



dp Z(z,p)Z(z,-p)[v + v ]x (p)], 



(A2) 



where Vi T (p) 







By using Z(z,p) = exp[— p 2 z/2] 



we perform integration only over an outer shell of radii A 
and A(l — 5), where 5 is a very small parameter. This is 
done essentially to integrate out the small-scale fluctua- 
tions. In the large- A limit, the contribution of the above 
term after the momentum-shell integration is 

A^5K d [v + v h (A)] 

By taking this term into account, one can express Z s 
viewed at a larger length scale (A — <5A) _1 in terms of Z s 
viewed at a smaller resolution as 

Z s (A - 5 A) = Z s (A) - Z s (A) K d A-< [v + v h (A)] , ( A4) 

where only the cutoff dependence is shown explicitly. 
This equation can further be transformed into a differen- 
tial equation form 



din Z s 
dA 



where L = A . By using the fixed- poin t value for u, one 
obtains the exponent given in Eq. (pM). 
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